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TURBULENT-HEAT -TRANSFER MEASUREMENTS 

AT A MACH NUMBER OF 2.06 

By Maurice J. Brevoort and Bernard Rashis 


SUMMARY 


TurLulent-heat-transfer measiirements were obtained through the use 
of an axially symmetric annular nozzle which consists of an irma-r shaped 
center body and an outer cylindrical sleeve. Measurements taken along 
the outer sleeve gave essentially flat-plate results that are free from 
wall interference and comer effects for a Mach number of 2.06 and for a 

Reynolds number range of I.7 X 10^ to 8.8 X lo'i^. The heat-transfer 
coefficients are in good agreement with available data from V-2 rockets 
and also check the Rubesin theory and the Van Driest theory for a Mach 
number of 2.0 and for a ratio of inner-surface to free-stream temperature 
of 1.8. The temperature -recovery factors are approximately O .5 percent 
lower than the factors given in NACA TN 2077 for a Mach number of 2.k. 



OOtbMbS 


INTRODUCTION 


The design of supersonic aircraft and missiles requires engineering 
information about heat-transfer coefficients and temperature -recovery 
factors for supersonic speeds that extend over a wide range of Reynolds 
niomber. In reference 1^ local -heat-transfer-coefficient measurements were 
presented for a Mach number of 5-05. Good agreement of these results with 
theoretical and experimental work was obtained. This method of testing 
and reduction of data is readily adaptable to obtaining accurate measvire- 
ments over an extended range of both Mach and Reynolds numbers • 

The purpose of this investigation is to extend the work initiated in 
reference 1 to a Mach number of 2.06. The same type of apparatus and 
method of reducing the data were used in this investigation as were employed 
in reference 1. The range of Reynolds number for which measurements were 

obtained is from I.7 x 10^ to 8.8 X lo”^. The results cover a tempera- 
ture difference of approximately 20° at 40 seconds after starting to 
approximately 5 ° 120 seconds after starting. The average value of the 

ratio of inner-surface temperature to free-stream temperature TwAco 
was 1.8. 
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SYMBOLS 


specific heat of sleeve material^ Btu/lb-°R 
specific heat of air at constant pressure^ Btu/lh-OR 


acceleration due to gravity^ ft/sec^ 

heat-transfer coefficient, Btu/ft^-sec-°R 

heat conductivity, Btu/ft-sec-'^R 

Mach number 

Nusselt number, hx/k 

Prandtl number, iiCpg^k 

Reynolds number, pVx/n 


Stanton number, ■ = — — — 

R Pr pVCpg 

average v ai l temperature, °R 

effective stream air temperature at vail, some tempera- 
ture which gives a theimal potential which is independent 
of heat-transfer coefficient h, 

stagnation temperature, °R 

inside -surface temperatirre of nozzle sleeve, °R 
free-stream temperature, °R 
time, sec 

free-stream velocity, ft/sec 

specific weight of sleeve material, Ib/sq. ft 

longitudinal distance along sleeve, ft (unless indicated 
otherwise ) 


X 



NACA TN 557^ 


5 


p 


T - T 

recovery factor, -S 2° 

T . - T 

t oo 

dynamic viscosity coefficient, Ib-sec^sq ft 
free-stream density of air, slugs/cu ft 

APPARATUS AND METHOD 


The apparatus consisted of an axia3JLy symmetric annular nozzle which 
was directly connected to the settling chamber of one of the cold-air 
blowdown jets of the Langley gas dynamics laboratory. The nozzle had a 
shaped wooden center body and two outer sleeves . The first sleeve was 
constructed from 8 -inch-diameter, extra heavy, seamless carbon-steel pipe 
and the second was constructed from l/l 6 -lnch stainless steel which was 
rolled into a cylinder and welded; both were surfaced-machined inside 
and outside to a wall thickness of 0.588 inch and 0.060 inch, respectively. 
A detailed drawing of the apparatus is shown in figure 1 which gives the 
location of the thermocouples and the static-pressure orifices. Several 
static -pres sure orifices and thermocouples were located 90 ° apart around 
the sleeve to calibrate the flow for axial symmetry. The coordinates of 
the center body are given in table I. The center body was designed by 
using three-dimensional characteristics. (See ref. 2.) Station locations, 
shown in figure 1 , refer to the distance in inches from the nozzle mini- 
mum station. 

Details of the installation of the thermocouples and static -pressure 
orifices are shown in figures 2 and 5 * The thermocouples are located 
0.060 inch from the inner surface of the sleeve. The wires are 
copper-constantan (O.OlO inch in diameter); conduction along the wire 
length is negligible. As indicated in figure 2, the thermocouples are 
in intimate contact with the steel wall so that thermal resistance at 
this junction is negligible. 

The Mach number distribution was measirred both lengthwise and around 
the sleeve. The orifices located around the sleeve were used to check the 
alinement of the center body within the sleeve. The results are shown in 
figure 4 and are accurate to iO.Ol. 

The temperature -recording equipment consisted of three synchronized 
high-speed self-balancing Brown ElectroniK strip chart recorders, having 
a total of 58 switches, two of which were connected to the settling- 
chamber thermocouples and the others were connected to the sleeve 
thermocouples . 
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A typical test consisted of operating the nozzle at the desired 
pressure. The stagnation temperature starts at essentially room tempera- 
ture and decreases as the piping is cooled, as illustrated in figure 5 
where stagnation temperature is plotted against time for a settling-chamber 
pressure of 106 Ib/sq. in. gage. The wall temperature also starts at essen- 
tially room temperature and tends to approach the eq_uilibrium temperat\ire 
which is approximately 25 ° F below stagnation temperature. This variation 
is shown in figure 6 where wall temperatirre is plotted against ti 3 re for 
stations 12 and l8 for a settling -chamber pressure of 106 Ib/sq. in. gage. 
Test runs were made for settling -chamber pressures of 1 , 58, I06, and 
185 Ib/sq in. gage. The test at a pressure of 1 Ib/sq in. gage was made 
with a sleeve O.060 inch thick, and the test setup was evacuated to an 
absolute pressure of 1.0 inch of mercury. Except for the first 20 seconds, 
the pressures were maintained constant for each test run. The recorders 
were calibrated immediately before and after each run and were found to 
be accurate to tl° F. 


In figure 7 are plotted, for various times during the test run, the 
values of wall temperature against longitudinal distance along the cylinder. 
These values were losed to determine the rate of change of the longitudinal 
2 

d T 

conduction k — along the cylinder. Test results were taken only for 

dx^ 


the length of the cylinder for which 


k 


^^Tav 

dx^ 


was zero. 


BEDUCTION OF DATA 


Short test runs in which time histories of the wall temperature are 
obtained give the essential data from which heat-transfer coefficients 
and temperature -recovery factors may be determined. The reduction of 
these data requires a method in which the stagnation temperature, and 
accordingly the equilibrium temperature, may vary in an arbitrary manner. 

The Stanton number is calculated frcan 

St = ( 1 ) 

pVCpg 

and the heat-transfer coefficient is calculated from 


h 


wc 


T - T^ 
w e 


(2) 
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By definition, the recovery factor is 


^r 



(3) 


The method of reducing the data is simply to select a recovery factor 
and then obtain from equation (5). Substitute this value of Tg into 

equation (2) and obtain values of h for different heat-flow rates. These 
values of h are then substituted into equation (l). The true values of 
Tgj and St ajre obtained when St is constant with time (for dif- 

ferent heat-flow rates ) . Figure 8 shows the values used in evaluating the 
Stanton number and recovery factor at station 12 for a settling-chamber 
pressure of 106 Ib/sq in. gage. 

The values of specific heat and specific weight of the sleeve mate- 
rial were taken from reference 3 . The values of viscosity and Prandtl num- 
ber ( 0 . 71 ) for air were also taken from reference 3 and were based upon 
the inside-surface temperattire of the sleeve. The value of T^ chosen 

was at 80 seconds after starting. (This arbitrary choice of temperature 
can be made since the viscosity coefficient is practically a constant 
factor of the heat -conduction coefficient.) 


FREX:iSION OF EQUTFNENT 


The basic data consist of time histories of temperature at each 
thermocouple location. The recorders were calibrated immediately before 
and after each r\m and were found to be accurate to ±1° F. The thermo- 
couple wires axe 30 -gage (0 .010-inch-diameter ) ccpper-constantan and of 
such small diameter that conduction along their length can be neglected. 
The wires are insiilated except at the jimctlon where intimate contact 
with the sleeve keeps any thermal resistance to a minimum. 

The sleeve is enclosed in a vacuum jacket to avoid free convection 
and, for the temperature range considered, radiation effects range up to 
2 percent. The product of the specific weight and the specific heat of 
the sleeve material as used in equation ( 2 ) is accurate to ±5 percent. 

The value of — used in equation ( 2 ) is that of the actual 
dt 

thermocouple reading, whereas it should be the value associated with the 
average teonperature in the thickness of the sleeve at each point. During 
the test runs the inner -surface temperatures are lower than the average 
temperatures by approximately 1° to 1.5° at the beginning of the runs 
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(iA) seconds after starting) and become the same at the end of the runs 
(I50 seconds). This deviation causes the derivative of inner-surface 
temperature with time to be less than the derivative of average tempera- 
ture with time by approximately 1 percent. 

The analysis of reference was used to determine the effect of the 
sleeve diameter upon the flow characteristics. This effect was computed 
to be less than one-hundredth of 1 percent. 


RESULTS AND DISCUSSION 


Figure 7 shows the variation of wall temperature with longitudinal 
distance along the cylinder for a settling-chamber pressure of 
106 Ib/sq, in. gage. Over the test range, the wall temperatures are con- 
stant along X. In the relation 

wc dTg^^/dt 


w and c are constant and dTg^^y'dt is constant because T^ is constant. 

Therefore, if there is to be a variation of h with Reynolds number or 
X, Tg must vary with x. The value of Tg obtained for the test at a 

settling -chamber pressure of I06 Ib/sq. in. gage, evaluated at 80 seconds 
after starting, actually decreases approximately 1.5° • 

Figure 9 shows the variation of local Nusselt number with local 
Reynolds number. The value of x used in evaluating these numbers was 
based upon an x that was considered to be zero at the nozzle minimum 
station. A single line faired through all the data points has a slope of 
approximately 0.8; but if lines are faired through the points for each 
settling -chamber pressure, the slopes are slightly higher. It is reason- 
able to assume, however, that the higher stagnation pressures would pro- 
duce earlier transition and the zero value of x would move upstream 
with increasing stagnation pressure. Hence, if it is assumed that data 
from each settling-chamber pressure should merge into a continuous line 
(rather than a series of steps), x may be adjusted, for zero x-locations 
(the effective beginning of the turbvilent boundary layer) which would 
make all the test runs coincident. This adjustment lias been made in 
figure 10 by using values of x eq_ual to zero at 1.^ inches downstream 
of the nozzle minimum station for a settling-chamber pressure of 
1 Ib/sq in. gage and 2.00, b.40, and 6.00 inches upstream of the minimum 
for the chamber press^Ires of 38, I06, and I85 Ib/sq in. gage, respectively. 
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The Nusselt numbers were found to vary from 1,690 to 44,950 for the 

Reynolds number range of I .7 X 10° to 8.8 x 10' (based on adjusted 
zero x-locations ) . For con^iarison, the ciirves based upon the analyses 
of Van Driest (ref. 5 ) sJ3.d Rubesin (ref. 6) are shown. In the Van Driest 
analysts T^y/Too was considered to be 1.8 (the average value of the test 

results). For comparison with these references, the data were ccaaputed 
by using free -stream temperature to determine the density and the velocity. 
The wall temperature was used to determine the viscosity and the Prandtl 
number. The experimental values from the data of V-2 rockets (ref. 7 ) 
are also determined in this manner. The Mach number range for the V-2 data 
in flgire 10 is frcm 1.95 "to 2 .O 7 . The agreement between the results of 
this investigation and those of references 5 , 6 , and 7 is good. 

Figure 11 shows the variation of local temperature -recovery factor 

with local Reynolds number. The variation is from O. 89 I at R = I .7 x 10^ 

to 0.870 at R = 8.8 X 10 . The results are ccmpared with a curve which 
represents the fairing of the data of reference 8 . The data of reference 8 
are for a Mach number of 2.4. The res^alts obtained in this investigation 
are approximately 0.5 percent lower. Also included for comparison are the 

curves for recovery factor eq.ual to pP/3 

and The wall tempera- 

ture was used to determine the Prandtl number. 


CONCLUDING REMARKS 


Turbulent -heat -transfer measurements that gave essentially flat -plate 
results were obtained for a Mach number of 2.06 and for a Reynolds number 

range of I .7 x 10^ to 8.8 x lo'^. The Nusselt numbers are in good 
agreement with theoretical analyses and with experimental data obtained 
with the V-2 rockets. Bie temperature-recovery factors are approximately 
0.5 percent lower than other available data for a I<fech number of 2 . 4 . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
langley Field., Va., December 5^ 195^* 
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TABIE I.- CENTER-BODY COORDINATES 


X, in. 


- 10.25 

-1^.7 

- 4.5 

-4.0 

- 3.0 

-2.5 

- 2.0 

-1.5 

- 1.0 

-.8 

-.6 


Radius, in. 


2.000 

2.000 

2.020 

2.100 

2.500 

2.735 

2.970 

3.150 

3.300 

5.340 

3.375 

3.400 

3.425 

3.4375 

5 . 456 A 

5.4335 

5.4280 

3 . 420 A 

3.4108 

3.3990 

3.3852 

3.5696 

3.3522 

5.5336 

3.3141 

5.2941 

3.2739 

3.2558 

3.2541 

3.2149 

3.1964 


X, in. 


Radius, in. 


5.1787 

3.1618 

3.1456 

3.1305 

3.1157 

3.1020 

3.0890 

3.0769 

3.0655 

3.0550 

5.0455 

5.0364 

5.0285 


2.8824 
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figure 1.-- Test arrangement. Dimensions are in inches. 
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Figure 2.- Thermocouple installation. 
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Figure it.- Mach number dletrlhutlon for settling- chamber pressure of 

106 Ib/sq in. gage. 
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Figure 6.- Variation of -wall temperature -with tiioe for a sett lin g" chamber 

flnTS=> nf* 1 Th/cn -^-n a»ocrO . 
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Figure 7.- Yariation of wall ten5>erature with longitudinal distance for 
a settling- chamber pressure of IO6 Ib/sq In. gage. 
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number for corrected locations of x = 0. 
numiber detemolned for wall ten^jeratures . 
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Reynolds number, R 

Figure 11.- Variation of local recoveiy factor with local Reynolds 
number for corrected locations of x = 0. Viscosity determined 
for wall teu^ieratureB . 
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